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Abstract. Recent developments of theoretical studies on spin nanotubes are 
£^ ■ reviewed, especially focusing on the S = 1/2 three-leg spin tube. In contrast to the 

three-leg spin ladder, the tube has a spin gap in case of the regular-triangle unit 
cell when the rung interaction is sufficiently large. The effective theory based on 
the Hubbard Hamiltonian indicates a quantum phase transition to the gapless spin 
liquid due to the lattice distortion to an isosceles triangle. This is also supported by 

£~0 ' the numerical diagonalization and the density matrix renormalization group analyses. 

^^ . Furthermore, combining analytical and numerical approaches, we reveal several novel 

magnetic-field-induced phenomena: Neel, dimer, chiral and/or inhomogeneous orders, 
new mechanism for the magnetization plateau formation, and others. The recently 
synthesized spin tube materials are also briefly introduced. 
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1. Introduction 

The geometrically frustrated low- dimensional quantum spin systems have attracted 
increasing attention in recent years. Among them the spin nanotube is one of the most 
interesting systems, because it is expected to be a new generation of multifunctional 
devices. Recently some candidates for the spin nanotubes have been synthesized; three- 
leg tubes [(CuCl 2 tachH) 3 Cl]Cl 2 [H|2] and CsCrF 4 [3], a nine-leg one Na 2 V 3 7 [4], and 
a four-leg one CU2CI4 • D8C4SO2 (Refs. [5] and [6]). Several theoretical works indicated 
various exotic quantum phenomena of the spin tubes [7H2TJ. Most of them concentrated 
on the S = 1/2 three-leg spin tube, because both of frustration and quantum fluctuation 
are the largest. In this paper, we review the recent theoretical results of the system. 

In the S = 1/2 three- leg antiferromagnetic spin tube the unit cell consists of 
three spins. According to the Lieb-Schultz-Mattis theorem [22J, the spin gap must be 
accompanied with at least doubly-degenerate ground states. In fact, previous numerical 
analyses [HIIIHKIT] have confirmed such doubly-degenerate S = ground states due to 
the spontaneous breaking of the translational symmetry along the leg direction. The 
ground states have a valence-bond type (dimerized) order [HHH]. Some numerical 
works [71[T6l[l7] indicated that the spin gap is quite fragile against the lattice distortion 
changing the unit cell from a regular triangle to an isosceles one. As a result, when one of 
the three rung couplings is changed, a quantum phase transition occurs from the spin gap 
to gapless phases. Developing the effective theory based on the Hubbard Hamiltonian, 
using the numerical exact diagonalization and the density matrix renormalization group 
(DMRG) approach, the ground-state phase diagram of the quantum phase transition 
was presented [7], reviewed in sections 3 and 4. 

This system was also theoretically revealed to exhibit some exotic phenomena in 
magnetic field. The numerical diagonalization study [12] suggested that a magnetization 
plateau appears at 1/3 of the saturation for sufficiently large rung interactions. 
Introducing the above lattice distortion, two different plateau formation mechanisms; 
the up-up-down (uud) and the dimer-monomer types, are expected to appear, depending 
on types of isosceles triangle. Our recent work [23] has indicated a new plateau phase 
with a dimer and/or chiral order accompanied with the staggered moment. On the 
other hand, the bosonization for the weak rung-coupling regime [18] has shown that 
a vector chiral phase is widely present in the field-induced Tomonaga-Luttinger liquid 
(TLL) phase. In addition, a spin-wave type approach p[9] has predicted that not only the 
chiral order but also an inhomogeneous magnetization order appears near the saturation. 
These field-induced phenomena will be discussed in detail in section 5. 

The valence bond solid (VBS) ground states were revealed to appear in some 
different spin tube systems. For example, the quantum Monte Carlo simulation study 
[24] suggested that several different types of the VBS ground state occur in some 
chiral spin nanotubes without frustration. The recent non-linear a model and DMRG 
approach [25] indicated some three-leg higher spin tubes exhibit several Haldane phases. 

One of the most realistic candidate of the S = 1/2 three-leg spin tube is 
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[(CuCl 2 tachH) 3 Cl]Cl 2 [2]. In this material each triangle unit cell is oriented upside- 
down to adjacent one, different from the straight spin tube (see section 6). The DMRG 
calculation of such a twisted spin tube JEHHIEE] suggested that the material has no spin 
gap, which is consistent with the magnetization measurement [2j. 

As a future prospect, we consider the superconductivity in the carrier-doped spin 
nanotubes in section 7. A possible chirality induced superconductivity will be proposed 
as a new mechanism. 

2. Isosceles spin tube 

We consider the 5 = 1/2 isosceles three- leg spin tube [7] shown in figure [Q described 
by the Hamiltonian 

3 L 2 L L 

7~L = h 2^ 2_^ S i,j ' S i,j+1 + ^ 2^ 2_^ Si ' j ' Si+hj + J ' r ^ S3 ' j ' Sl,j ^ 
i=l j=l i=l j=l j=l 

where Sij is the spin- 1/2 operator and L is the length of the tube along the leg direction. 
The exchange interaction constant J\ stands for the neighboring spin pairs along the 
legs, while J r and J[ the rung interaction constants. All the exchange interactions are 
supposed to be antiferromagnetic (namely, positive). The ratio a = J'J J r expresses the 
degree of the asymmetry of the rung interactions. Throughout this paper, we fix J T to 
unity. The effect of the magnetic field H along the z direction is taken into account by 
adding the Zeeman term "Hz = —HJ2ij $i,j to (00). 



Figure 1. Sketch of the isosceles three-leg spin tube. 



3. Effective theory 

In this section, we explain low-energy effective theories for the spin tube ([I]). In 
section I3.1[ we sketch a method to draw a phase diagram in the whole coupling-constant 
space (a,Ji) used in [7]. Next, we investigate two special regimes J r ^C J' t ^ J\ and 
J r 3> J\, respectively, in sections 3.1 and 3.2. 

3.1. Global phase diagram derived from the Hubbard model on the tube lattice 

Here, we explain a systematic method to draw global phase diagrams of one-dimensional 
antiferromagnetic quantum spin systems. It is well-known that S = 1/2 Heisenberg 
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model on an arbitrary lattice is obtained from the corresponding half-filled Hubbard 
model in the limit of strong on-site Coulomb interactions. Especially in one dimension, 
the spin configurations of the low-energy states in the Heisenberg model agree with 
those in the half-filled Hubbard model even with a weak Coulomb interaction [261428] . 
The phase in the weak- Coulomb regime often smoothly connect to those in the strong- 
Coulomb regime in one-dimensional electron systems. By these arguments, we obtain 
a low-energy effective theory for the spin tube ([1]) from the corresponding Hubbard 
model. To discuss the wider parameter space, first we diagonalize the kinetic parts of 
the Hubbard Hamiltonian including both the leg and rung hopping terms [29H3T] . Then, 
we take account of the on-site Coulomb interaction as the perturbation, with help of 
the Non-Abelian bosonization [32H34"] and conformal field theory (CFT). 
The Hamiltonian of the Hubbard model on the three-leg tube lattice 

n. — %hop + 'Hint (2) 

consists of the hopping part 

L 3 
rihop — / J / J / J {^ C n+l,i,a C n,i,(T + ■Si+l.iC ?lj j +lj0 .C n) j )0 - + II. C.J 
n=l j=l <T=t,4- 

and the on-site interaction part 

L 3 

%int = U E E U ^M n n,i,l (3) 

n=l i=l 

where n n ^^ = c* nia c n ^ a and U > is the repulsive coupling constant. The electron 
operators c n> i >a and c\ { a satisfy the periodic boundary conditions for both the leg and 
the rung directions, c n+ i,i,a = c n ^ a , c nt i + ^ a = c n ^ yCT) and anticommutation relations. 
The hopping parameters are given by t > 0, si,2 = £2,3 = s > and 33,1 = (3s > 0. The 
strong coupling expansion shows that this model at the half-filling case is reduced to 
the Heisenberg model with J x = At 2 /U, J v = As 2 /U and a = /3 2 . 

By a suitable unitary transformation, the hopping Hamiltonian can be written into 
the following diagonal form 

3 

^hop = E E E E i( k ) d ii,*dk,i,e (4) 

k i=l o-=t,J. 

where the wave number k is summed over ^ < k < 2tt. With an orthogonal matrix 
Oij, the operators dk,i, a and d kicr are defined by 

1 l 3 

d kM = —= E E eX V(- ikn )°ij C n,j,a, (5) 

V L n=l j=1 

1 L 3 

d U,a = 7>= E E e M^n)O ljC i tJ(T (6) 

n=l j=l 
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which satisfy the standard anticommutation relations. The energy eigenvalues of the 
one-electron states are 

E 1 (k) = -ps + 2tcoak, (7) 

Mk) = \{ps - s^/3 2 + 8 + 4tcos k), (8) 

E 3 (k) = -(/3s + s^/3 2 + 8 + Atcosk). (9) 

Note that a degeneracy E\(k) = E 2 (k) appears at (3 = 1 due to the permutation 
symmetry (S 2 ,j -H- Sij) for i — 1, 2, 3. 

For the half-filled case, 3L one-electron states should be occupied by electrons with 
up and down spins. As a result, the ground state of the hopping Hamiltonian has one, 
two or three pairs of the Fermi points (kj,kj) just on the Fermi sea, depending on 
the parameters s/t and f3. Since the low-energy excitations are given by the particle- 
hole creations around these Fermi points, they can be represented by using the Dirac 
fermions, the left mover ipj jCr (x) and the right one ipj^(x), which are defined from the 
electrons around the j-th pair of Fermi points. If the j-th band has no Fermi points in 
the half-filled case, we should neglect ipj ttr (x) and ^^(x). On this understanding, we 
approximate the original electron operators in terms of the Dirac fermions as follows: 

3 

c n ,i,a ~ V^^ °^ ( ex P(% a; / a )'0j>(a;) + exp(ifc j x/a)^ J> (a;)} (10) 

i=i 
where a is the lattice spacing with dimension of length and x = an is the continuous 
position coordinate. We add the on-site Coulomb interaction (|3]) to this free Dirac 
fermion system, and use the Non-Abelian bosonization techniques. Following the field- 
theory argument in [28], we expect that for the number of Fermi-point pairs is odd 
(even), the spin excitations are gapless (gapped) in the half-filled Hubbard tube. In 
particular, in the cases of one or two Fermi-point pairs, we can explicitly determine 
whether or not a spin gap exists as follows. 

First, we consider the case of one pair of Fermi points k 1 = ^ and k\ = |. In this 
case, the interaction ([3]) is approximated as the sum of a Umklapp interaction and two 
marginal ones 

Hut ~ dx flfiOi (ar) - g%0 2 {x) - g^O z (x) + • • 

where #1,2,3 are positive coupling constants proportional to U. The Umklapp term is 
expressed as 

1 {x) = ^ lit (x) t ^i,;(x) t ^i, t (x)^i,i(x) (12) 

and the marginal interaction between the U(l) charge currents is given by 

2 (x) = MrfMx)iJ>i(.xNi(?) (13) 

where ipi = t (tpi^,ipi t i). It is known that the bosonized form of 0\^ 2 contain only the 
charge degrees of freedom and they open a charge gap when g 2 is positive. Then, 
the remaining spin degrees of freedom are described by the gapless level-1 SU(2) 



:n} 
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Wess-Zumino-Witten (WZW) theory [32-34J. This phenomenon, i.e., the charge-spin 
separation is well-known in the single Hubbard chain model. For this WZW theory, the 
third interaction 



<D 3 (x) = i)i{x) ] — ^i(ar) ■ i)\{x) ] '—ipi(x) 



(14) 



is known to be marginally irrelevant if g^ > 0. Except the above interactions 0\^^{x), 
there is no relevant operator with the invariance under the one-site translation along 
the leg, 



VvO) ->■ e faA x )> ^vO) -► e ^iA x ) 



(15) 



as in the case of the Heisenberg chain. The spin excitations remain gapless. 

On the other hand, when there exist two pairs of the Fermi points, (ki,ki) and 
{k2,k 2 ), the spin excitations are suffered from relevant interactions. In this case, the 
spin sector in the hopping part of the Hubbard tube is described by a level-2 SU(2) 
WZW model derived from two decoupled Dirac fermions [38J. The Coulomb interaction 
yields several perturbations for this theory. For example, the following term 

^i,t(a ? ) t ^2,t(a ? )^u(»)V2 I 4.(a?) ( 16 ) 

contains a relevant perturbation in the level-2 WZW model, which is invariant under 
the translation 



il)jA x ) -> e j ^, CT (x), il)jA x ) ~> e J ^jA x ) 



(17) 



Particularly for — 1, much more relevant operators are allowed by the translational 
invariance because of the coincident Fermi points k\ = k 2 and k\ = k 2 . Therefore, 
we conclude that any gapless spin excitation generally has no chance to survive except 
particularly rare cases [e.g., when the coupling constant of (fT6|) is zero]. For the case of 




Figure 2. Phase diagram obtained from the effective Hubbard model @. In the 
strong-coupling limit (U/t,U/s 3> 1), the horizontal axis f3 2 and vertical one (s/t) 2 
can be regarded as a and J r /Ji, respectively. The effective theory claims that the 
phases (II) and (IV) are gapless and the phase (III) is gapful. The phase (I) is possibly 
gapless. 
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three Fermi-point pairs, the interactions among three Dirac fermions are generated from 
(J3J), and thus it is difficult to analyze them and judge whether or not the spin excitation 
can survive as gapless. The gapless spin excitation is expected from [281135] . 

From these arguments, we can draw the ground-state phase diagram of the half- 
filled Hubbard tube shown in figure El The phase (I) has three pairs of the Fermi 
points, the phases (II) and (IV) have one pair, and the central phase (III) has two 
pairs. Therefore, we predict that the phases (II) and (IV) have gapless spin excitations, 
whereas the phase (III) possesses a spin gap. By (t/s) 2 = J x j J\ and f3 2 = a obtained 
from strong on-site interaction, figure [2] depicts the phase diagram of the S — 1/2 three- 
leg spin tube (CQ). Since we have treated the on-site interaction perturbatively in the 
effective theory, we should not trust the obtained phase boundaries as accurate ones. 
Note that the gapful phase (III) is predicted to be extended around the line (3 = 1 for a 
finite s/t. In the limit s/t — > oo, the both the left- and right-side phase boundaries of 
the region (III) converge to (3 = 1. This narrowing of the phase (III) is consistent with 
the numerical results (HUE] in the strong-rung-coupling limit J T /J\ — > oo. 

Finally, we argue the universality classes of the phase transitions at two phase 
boundaries, (II)- (III) and (III)- (IV). For the level-1 SU(2) WZW model in the phases 
(II) and (IV), the most relevant perturbation is the marginal current- current interaction 
[32Tf3l] . (Tbf|) . which is the only allowed invariant interaction. Since it can be marginally 
"relevant" when parameters are finely tuned and then g% becomes negative, we speculate 
that the transition from the phases (II) or (IV) to (III) is caused by this marginal term. 
Therefore, the transitions are expected to be in the Berezinskii-Kosterlitz-Thouless 
(BKT) universality class 



3.2. Gapless phase for J r <C J' v <C J x 

In section 13.11 we have obtained a qualitative phase diagram of the spin tube (CQ) as 
shown in figure [2j However, some subtle points still remains. Particularly the existence 
of the gapless phase (I) should be discussed. Here, we focus on the right-lower regime 
in figure [2] J r <C J' x <C J 1; [a = J' r /J r 3> 1). In this regime, we discuss the existence 
of the gapless phase (I) and the phase transition between the gapless phase (I) and the 
gapful phase (III). 

We introduce three level-1 SU(2) WZW models for the three decoupled Heisenberg 
chains [32TI31] . and we treat their rung couplings as the perturbation The first and the 
third chains are coupled to each other with J[ which is much stronger than two remaining 
couplings J r . It is well-known that the J' r coupling involves a relevant interaction with 
conformal dimensions (1/2, 1/2) in the two coupled WZW models, it produces an energy 
gap [39J. This relevant operator behaves like the mass operator in the corresponding 
effective field theory. On the other hand, like the case of one Fermi-point pair in section 
13. H the WZW model for the second chain has the marginal irrelevant interaction — g^O^ 
with a finite negative coupling constant —g% < 0. The operator O3 is equivalent to ( 1T4l) . 
if we use the Dirac fermions (^1,0-, ^l <x) to describe the second chain. The negative 



Quantum spin nanotubes — frustration, competing orders and criticalities - 8 

sign makes —g^O^ irrelevant and the second chain is gapless. A weak rung coupling 
J r between this WZW model and the massive theory for the two coupled chains must 
give a correction to the coupling constant —g^. If J r is sufficiently small, the sign of 
— g 3 < would not change and the gapless excitation is preserved. These arguments 
convince us that the gapless phase (I) definitely exists. Furthermore, we expect the 
phase transition from the gapless phase (I) to the gapful phase (III). Namely, if J r 
exceeds a critical value, the coupling constant —g 3 might change to be positive. In this 
case, the operator — g-^O^ becomes marginally relevant, which produces an excitation 
gap. Since the operator is marginal, this transition between the phases (I) and (III) is 
of the BKT type [36,37J, as in the zigzag Heisenberg chain [ID]. The effective theory 
for J T <C J' T -C J\i (a = J' r /J T ~^> 1) is obtained, and the transition from the phase (I) 
to (III) is discussed in [7]. 

3. 3. Perturbation theory from the strong rung coupling limit 

In case of J r 3> J\, the perturbation theory from the independent triangle limit will 
work well. Let us begin with the J\ — case (i.e., three-spin problem), which can 
be easily solved and there are 8 eingenstates. Here ipQ (S to t, S^ ot ) represents the wave 
function with eigenvalues S to t, Sf ot an d P, where P is the eigenvalue of the rung-parity 
operation Sij <^> S^j. The state | til2T3) is abbreviated as | tlT) f° r example and [i,j] 
denotes the singlet pair [i,j] = (1/a/2)(| tilj) ~ I Utj))- We show four eigenstates with 
S^ ot > only, because other four eigenstates with S* ot < can be easily obtained by 
interchanging t^4-j 



/ o o \ 


state 


energy 


'-'tot 


qz 
°tot 


P 


* +> (§- + §)= 


: | ttt) 


1/2 + a/4 


3/2 


+3/2 


+1 


#ii4 = 


= ^(im> + itit> + utt>) 


1/2 + a/4 


3/2 


+ 1/2 


+1 


« +, G-4)- 


: i=(im)-2int) + iitt)) 


-l + a/4 


1/2 


+ 1/2 


+ 1 


■^{IA)- 


= -^(im)-utt» = it> 2 [i,3] 


-3a/4 


1/2 


+1/2 


-1 



;is) 



These wave functions do not depend on a because all of them are completely classified 
by the eigenvalues S tot , S* ot and P. The lowest energy states are ipQ (1/2, ±1/2) for 



a < 1, while they are ipQ \l/2, ±1/2) for a > 1. 



The first order perturbation theory with respect to Ji, retaining four states with 
Stot = 1/2 having lower energies and neglecting other 4 states with S to t = 3/2 leads to 

U cS = ^Y, T rT J+ i (l + 8[^ +1 + a^ +1 ]) - (1 -a)$>? (19) 
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where T is a spin-1/2 operator acting on S* ot and a is also a spin-1/2 operator acting 
as 

* + 4 + \sz*) = o, cr+4-\s? ot ) = 4 +) (^ ot ) (20) 

*-4 + \s? ot ) = 4~\s? ot ), *-4~\sz ot ) = o 

where S tot = 1/2 indices were omitted for simplicity. The energy difference between 
Vt) Wtot) an d V'o Wtot) ac ^ s as a "magnetic field" applied on a spins. This W, e g is 
essentially the same as that obtained by Nishimoto and Arikawa [T7| . 

In case of sufficiently strong asymmetry \a — 1| ^> 1, the a spins are completely 
polarized and only T spins survive, where 7i e s is reduced to that of the antiferromagnetic 
Heisenberg chain having the gapless excitation. Thus the phase (III) has a finite width 
along the line J r / J\ = const ^> 1. Next, let us focus on the region around the symmetric 
line a ~ 1. In the limit J\j J t — > 0, the ground state has saturated cr spins at any 
asymmetric point a/1. For a finite Ji/J T , however, we expect an extended gapful 
phase around a = 1, if the energy gap exists at a = 1 due to the coupling between T 
and cr spins. The finite energy gap does not vanish by an infinitesimal external field 
a — 1. In other words, the magnetization process of the o spins should show a zero- 
magnetization plateau. Therefore, an energy gap would also be present for sufficiently 
weak asymmetry \a — 1 1 ^C 1 . 

For the rest of this subsection we suppose a = 1 (i.e. regular triangle case) where 
the translational invariance by one site along the rung direction holds and the lowest 
energy states are 4-fold degenerate. In this case another set of useful expressions for the 
lowest energy states is 

1 1 L) = -^ (| m> + w| tit) + uT 1 ! ttt» (21a) 

| t R) = -^ (| ttl> + ^1 tit) + w| |tt)) (2U) 

1 1 L) = -L (| ||t) + w| It!) + u~ l \ t!!)) (21c) 

\iR) = l= (| ||t) + w -i| iU ) + oj\ HI)) (21d) 

where a; = exp(27ri/3) and indices L and i? denote the wave number A; = 2tt/3 and 
fe = — 2tt/3 along the rung direction, respectively. The first order perturbation theory 
with respect to Ji, retaining above four states with S to t = 1/2 and neglecting other 4 
states with S to t — 3/2 leads to [TOlITT] 



n *s = f E T ; • T ^+i C 1 + A ^t r I + i + T^i+i)] ( 22 ) 

j 

where T is a spin-1/2 operator acting on the first indices of the above states and t^ 1 are 
the spin-1/2 matrices acting on the second indices as 

r+|L) = 0, r-\R) = \L) 

t+\R) = \L), r-\L) = [A6) 
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Of course this "H e ff can be obtained from (1191) through the unitary transformation and 
letting a — 1. 

Schulz [in] analyzed "H e ff using the Jordan- Wigner transformation and concluded 
that the system is gapped. Kawano and Takahashi [11] performed the DMRG calculation 
for Hcfi to find that Hamiltonian 7i e s has a spin gap of 0.277Ji. When the next- 
nearest- coupling term having a special coupling constant is added to fl22l) they found the 
exact ground state closely similar to the Majumdar-Ghosh state, which is spontaneously 
dimerized and breaks the translational symmetry. They insisted that the ground state of 
( j22|) is similar to this Majumdar-Ghosh state, and provided strong numerical evidences. 

4. Spin gap and ground-state phase diagram — numerical study 

In this section the numerical ground state phase diagram of the isosceles spin tube 
obtained in our previous work [7J is presented. The method to derive it is the 
phenomenological renormalization combined with the numerical diagonalization up to 
L = 10 and the DMRG up to L = 128. 

A useful order parameter to determine the phase boundaries between the gapless 
and the gapful phases is the spin gap A, which is the energy gap between the singlet 
ground state and the triplet excited state in the finite, but large system. We calculate 
it by means of the DMRG up to L = 128. 

In order to determine a phase boundary of the usual second order phase transition, 
the phenomenological renormalization equation LiA Ll (a c ) = L 2 A L2 (a c ) is often used 
effectively. For the present critical point a c , however, this type of phenomenological 
renormalization has no clear crossing point, but the scaled gap LA increases with 
increasing L in both gapless and gapful phases. This is because the scaled gap is 
an increasing function with respect to L not only in the gapful phase but also in the 
gapless phase, since the finite-size gap must have the logarithmic size correction term 
~ — 1/logL. Here, we should recall that the logarithmic correction normally vanishes 
just at a c due to the SU(2) symmetry in the c = 1 CFT |40j. Therefore, instead of 
using the crossing point of the scaled gaps, we can estimate a c as a point where the 
size correction is minimized. Note that the minimum value L 1 A^ 1 — L 2 A_l 2 decreases 
as the size increases. This phenomenon and the assumption of the BKT type transition 
suggest that this minimal value approaches zero as the system size increases. This is 
quite reasonable if we suppose the most important finite-size correction to the scaled 
gap LA next to 1/logL term is order of 1/L 2 [4*0144*2] . We thus determine a c from the 
minimums for two large systems with L\ = 96 and L 2 = 128 for J x < 2. The estimated 
a c \ and a c2 are shown as crosses in figure El They correspond to the phase boundary 
between two regions (II) and (III) and that between (III) and (IV), respectively. At 
least these boundaries for the strong-rung-coupling regime J\ < 2 are precise enough 
to justify that a finite gapful phase (III) exists. However, it is difficult to obtain a c for 
J\ > 2, because the DMRG calculation does not well converged there. 

In order to determine the phase boundaries for the weak-rung-coupling regime 
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Figure 3. Ground-state phase diagram of the isosceles three-leg spin tube (Q]), derived 
from the numerical analysis. The phases (I)-(IV) correspond to those in figure [2] 
The cross points are determined by the DMRG, and the solid lines by the numerical 
diagonalization. The circle points are obtained from the level-spectroscopy method in 
section 4. 



J\ > 2, we use the minimum points of LiAl ± — L 2 Al 2 calculated by the numerical 
diagonalization up to L = 10 under the periodic boundary condition. Using the 
estimated phase boundaries for (Li,L 2 ) = (6,8) and (8,10), and assuming the size 
correction is proportional to 1/L 2 in both directions of J\ and a the phase boundaries 
among the phases (I), (II), (III) and (IV) in the thermodynamic limit were estimated. 
The phase boundaries are also shown as solid curves in figure El At least the phase 
boundaries (II)-(III) and (III)-(IV) are consistent with the DMRG results for J\ <C 1. 
The boundary (III)- (IV) is, however, significantly deviated from the DMRG estimation 
for J\ ~ 1. This discrepancy is supposed to be due to the error of extrapolation. This 
analysis also justifies the existence of the phase (I). However, the error of extrapolation 
becomes larger as we approach the line 1/ J\ = in the case of a < 1. Thus it is difficult 
to conclude that the phase (I) really exists for a < 1 within the present numerical 
demonstration. It is also difficult to confirm the boundary (I)- (II), and the phase (I) 
might combine with the phase (II) in a certain regime with a < 1. 

The phase boundaries between the phases (I) and (III), as well as (III) and 
(II), for smaller Ji (0 < J\ < 0.5), were well reproduced by the level spectroscopy 
method [10~ tH3lfi5] , which is one of precise techniques to estimate the critical point for 
the BKT transition. But the level spectroscopy method does not work well for larger 
Ji, because of too large finite-size corrections. 

In our previous work [7], applying the conformal field theory analysis [16H51] . the 
central charge and the critical exponents of the spin correlation functions were estimated. 
The results suggested that all the phase boundaries in figure [3] belong to the BKT 
universality class. 
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5. Field-induced phenomena 

In this section, we consider effects of applied magnetic field in three-leg spin tubes. 
Particularly, we focus on field-induced phases (vector chiral order and magnetization 
plateau) and quantum phase transitions. 

5.1. Vector chiral phase in the weak rung-coupling regime 

Here, we study the weak rung-coupling regime in a magnetic field [18], namely, the 
Hamiltonian, 

H H = U + Hz (24) 

with Jx ^> |J r |, \J' T \. If a sufficiently strong magnetic field is applied and a finite 
magnetization occurs, the spin-rotational SU(2) symmetry is reduced to the U(l) type, 
in which the Abelian bosonization [33 | l34" | l52] becomes reliable at least for the weak 
rung-coupling regime. At the zero rung-coupling limit, the low-energy properties of i-th 
spin chain in the tube can be described by a free boson theory whose Hamiltonian is 
given by 



H, 



cff 



Jdx^K-^hy + KidJi) 2 ] (25) 

Here, <f>i(x) and 0i(x) are the pair of dual scalar fields (x = ja ), and K and v 
respectively denote the TLL parameter and the excitation velocity. The spin operator 
is also bosonized as 

S z hj « ^=d x <p t (x) + (-l) j A cos(v / 4T^ + 2ttMj) + ■■■ 

\/7l 



S£. « exp(iV^6i)[(-iy B + B x cos(V47r0 4 + 2ixMj) + •••], (26) 

with A n and B n being nonuniversal constants. Here, M = (Sfj) is the magnetization 
per site. In the present notation, K runs from 1/2 to 1 when the magnetization 
M is increased from to the saturated value 1/2. Using the formula ( 126]) . we can 
obtain the bosonized form of the perturbative rung coupling, and the resultant effective 
Hamiltonian for the spin tube is expressed as 



H cS = I dx^ V -\K-\d x $ q ) 2 + K{d x ® q ) 

J q=0 

+ 2 x -MJ r d x % + -^={J r - j;)(2^.$ + V2d x $ 2 ) 

V 7T A/37T 

+ ^( 2 j r + j;)(9,$o) 2 -^j;(9,$ 1 ) 2 

- ^(2j r - 1)(^.$ 2 ) 2 - y^(j v - j' T )d x %d x <$> 2 

+ Bl j r0 Q V(e 1; e 2 ) - b 2 (j t - jfao 1 cos(v / 27e 1 ) 

+ Y •W^(2$i, 2$ 2 ) - -y ( J r " J rK* cos(2v / 2T$i) + • • • . (27) 
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Here we have neglected all the terms with oscillating factors exp^Mnj) or (— 1) J ", 
and have introduced new pairs of boson fields; ($ , ©o) = (Si=i &■> Si=i 0i)/\/3, 
($i,6i) = (0i - 3 ,#i - 3 )/V2, and ($ 2 , @ 2 ) = (0! + 3 - 20 2 ,#i + # 3 - 20 2 )/VH. 
The first line in (127)1 is equivalent to three copies of TLLs in decoupled chains. The 
second line is linear terms which induce just a small correction to the magnetization M. 
In the isosceles case J r ^ J' r , (SfJ) = (S3 »•) 7^ (Sfj) would be realized. The third and 
fourth lines are quadratic terms changing the values of K and v. In the regular-triangle 
case, the quadratic part is diagonalized in the basis of ($0,1,2, ©0,1,2), and resultant TLL 
parameters and velocities for ($0,1,2, ©0,1,2) sectors, -^0,1,2 and f 0,1,2, are calculated as 
K = K{1 + ^a)" 1 / 2 , K l = K 2 = K{1- ^Jf 1 )" 1 / 2 = ' K g , v = v(l + ^f^) 1 ' 2 and 
V\ = v 2 = v(l — KJi0 -o y/ 2 = xig. The fifth and sixth lines correspond to the vertex- type 
perturbations, and the function V(e±, e 2 ) is defined by 



V(e h e 2 ) = 2 cos (^J2t^\ cos U/W/2e 2 \ + cos U/^Ke^j (28) 



In our notation, the scaling dimension of vertex operators are given as [exp(inv47T0j)] = 
n 2 K and [exp(iri\/7r0j)] = n 2 /(4K) at the decoupled case. We should emphasize that 
( 1271) does not include any relevant terms with $ and ©o except for commensurate cases 
with M = q/p (q andp: integer). This property is protected by U(l) spin-rotational and 
translational symmetries of the spin tube. Therefore, the ($0, ©0) sector is described 
by a TLL. On the other hand, the remaining sectors are subject to the vertex terms. 

We first focus on the regular-triangle case J r = J' r analyzing the above effective 
Hamiltonian ( 127J) . In this case, two vertex perturbations cos(v / 2~7r@i) and cos(2\/27r<I?i) 
in (127)1 vanish, and the scaling dimensions of V A (@i,@ 2 ) and V A (2$ 1 ,2$ 2 ) are readily 
evaluated: [V(Q 1 ,e 2 )] = U( 2K g) and [V(2$x, 2$ 2 )] = 2K g , In the weak rung- 
coupling regime, K g is very close to K and K > 1/2 holds in the magnetization 
process. Therefore, ^(@i,0 2 ) is more relevant and ($i, 2 ,0i, 2 ) sectors have gapped 
spectra. The form of the potential V^(Oi,0 2 ) is shown in figure H](a), in which the 
diamond region is meaningful in the full 0i-0 2 plane and it is called here physically 
relevant zone. Remarkably there are two minimum points in the diamond zone: 
(0i, 2 ) = (±\/27r/3, a/27t/3) = X±. A minimum is mapped to the other via sign 
change ©i — > — ©i that can be realized by the exchange of two chains S\j -H- S^j and 
Q\ -H- 6* 3 . This suggests a spontaneous breakdown of rung-parity symmetry. Vector 
chiralities in rung bond «y = (Sy x Si + ij) a is a natural candidate for the order 
parameter, and the bosonized forms of their z components are given as 



n z ld « - B 2 sin (v^T^i + v / W2©2) + ■ ■ 
n z 2tj « - B 2 sin (v 7 ^©! " \/3^7202) + • • 



k z 2J « B 2 sin ( V27T01 J + • • ■ . (29) 

At the point X + (A_), k\, becomes positive (negative). From the translational 
symmetry along the rung, we can predict (ftf »-) = ( K 2j ) = ( K 3j)- We thus conclude that 
in the weak rung-coupling regime with J r = J[ and a finite M, the low-energy physics is 
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V(0 0) Physically relevant zone \®l 




-2 -1 

Minimum points 



V a (0 1/ x/27c/T) «=> 5 




Figure 4. (a) (color online) Contour lines of potential F(@i,02) at the regular- 
triangle case, and (b) potential V a (@i, y27r/3) in the isosceles case with several values 
of a. 



governed by the TLL in the ($o? ®o) sector and the remaining gapped sectors generate 
a vector chiral long-range order with spontaneously breaking the rung-parity symmetry. 
In the commensurate case of M = 1/3, there appears the magnetization plateau at least 
in a strong rung-coupling regime (see the next subsection). It has been predicted in [12] 
that the plateau survives up to a fairly weak rung-coupling regime (J r ~ O.lJi) at the 
regular-triangle case J r = J[. From the bosonization viewpoint, the plateau is attributed 
to the emergence of a spin gap of ($o, ®o) sector induced by a higher-order perturbation 
cos(vl27r$o( a; ) + 6vrMj). In fact, this vertex term may be relevant at M — 1/3 in which 
the factor exp(±i67rMj) is equal to unity. This corresponds to so-called gapped chiral 
phase; the system has a gapped spectrum and a vector chirality is long-range ordered. 
In other words, the spin tube offers a unique chance to observe both gapless and gapped 
vector-chiral phases in the magnetization process. 

Next we consider effects of the rung deformation, J T ^ J' T . Even in this isosceles 
case, the chirality between the first and third chains (k 3 j) is still useful as a order 
parameter because the tube is invariant under rung parity operation between these two 
chains. Due to the rung deformation, all kinds of the rung coupling terms in ( 1271) are 
modified. As already mentioned, the boson linear terms just slightly change the uniform 
magnetization M. The quadratic term d x $>od x $>2 cannot be diagonalized in the present 
basis, but its effects would be quite small and the qualitative nature of the vector chiral 
phase is expected not to be changed. The most important point of the rung deformation 
is the change of the vertex potential due to an additional term cos(v27r@i) in the fifth 
line of (T2T|) . The modified potential is expressed as follows: 



v a (e 1 ,e 2 ) = v(G 1 ,e 2 



Jr 



COS 



27t6i 



2 cos (^V/2Q^\ cos (y37r/2e 2 ) + a cos (\/2^Q 1 



(30) 



Figure IHb) presents this potential as a function of 0i with fixing 2 = a/27t/3. It clearly 
shows that the potential form is changed from a double-well type to a single-well type 
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for a < 1, while the form is always double- well type for a > 1. This classical argument 
predicts the critical point a c = 0.5, and the chiral order is expected to disappear in 
a < a c . To gain deeper understanding of this expectation, we consider the effective 
theory for ($i,Oi) sector. If we naively replace cos(y37r/2@2) with its mean value 
C*2 = (cos(a/37t/262)) in the potential V a (@i, 62), the effective Hamiltonian is obtained 
as 

+ BlJ r a^ 2C 2 cos( v /^72ei) + acos(v / 2^6i) +■■■. (31) 

This is nothing but a double-frequency sine-Gordon model. Under the condition 
K\ > 1/4, this model is believed to exhibit an Ising-type quantum phase transition 
by tuning the ratio of two coupling constants aj^lC-i) [53j|54]- We can hence predict 
the phase diagram of the isosceles spin tube with a weak rung coupling and a finite 
magnetization M ^ 1/3 as shown in figure [5j The chiral ordering occurs in the region 



TL liquid + vector chiral order Q£ TL liquid 
►•« 



c 
,„v 1 rv 



I 



11 . .. i/a d 1 toe 

2-leg ladder Regular 

+ chain triangle 



Ising transition 3-leg ladder 



Figure 5. (color online) Ground-state phase diagram of the isosceles spin tube in a 
magnetic field in weak rung-coupling regime J\ 3> | J r | , \J' r \- 

< a < 1. The true critical value a c must deviate from its classical one 0.5. The 
order parameter (k^j) behaves as ~ (a — a c ) 1//8 near the transition point. From this 
figure, one also finds that once a coupling between a two-leg ladder and a single chain 
is introduced, a vector chiral order immediately emerges. 

Similar scenarios of a vector chiral order and an Ising transition are also expected in 
other tubes consisting of spin chains with a TLL parameter K > 1/2. For instance, even 
in the zero-field case, spin tubes with easy-plane XXZ anisotropy must yield a chiral 
order in the weak rung-coupling regime. It is known that the field-induced TLL phase 
satisfies K > 1 in a spin-1 antiferro magnetic Heisenberg chain [55] and the leading term 
of bosonized spin operator Sj~ has the same form as that of (1261) [561EZ]- Therefore, 
three-leg spin-1 tubes have a vector chiral order at least in high-field and weak rung- 
coupling regime. In the vicinity of the saturation field, three-leg spin-S 1 tubes can 
be analyzed by spin- wave type approach [19J. It also predicts the emergence of the 
vector chiral order in a certain weak rung-coupling regime. These results obviously 
suggest that a vector chiral order is generally induced by applied magnetic field in 
three-leg spin-S 1 antiferromagnetic tubes with a weak rung coupling. The spin-wave 
approach also points out the possibility of an inhomogeneous magnetization phase with 
(Sfj) = (^3,j) ¥" (^2j) [IS] m the regular-triangle case. 
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5.2. 1/3 plateau 

The magnetization plateau at M = M s /3 surely exists for J\j J T — )■ since the triangle 
unit cell is composed of three S = 1/2 spins, while it does not exist for the J\/ J T — > oo 
because the system is reduced to three independent S = 1/2 spin chains in this limit. 
Thus the M s /3 plateau phase diagram on the J\ — a plane is of immense interest. 

As discussed in section 13.11 we study magnetization plateaux in the effective 
Hubbard model under external field. The Zeeman term for electrons is given by 

L 3 



n z = -H/2j2Y,( c L 



t C »M>t C n,i,l 



Cn,i,i.) ) 



(32) 



which should be added to the Hubbard Hamiltonian. Elaborated analysis of the effective 
theory for this Hubbard model gives the conditions for the existence of the spin excitation 
gap. Then, we show the existence of the 1/3 magnetization plateau and nonexistence of 
any other nonzero magnetization plateaux. Details will published elsewhere [58]. The 
phase diagram of the 1/3 magnetization plateau is depicted in figure El This agrees with 
the numerical calculation qualitatively, as shown later. 

2.0 



a 



1.0 



Plateau 




Figure 6. The phase diagram of the 1/3-magnetization plateau predicted by the 
effective Hubbard model in the strong-coupling limit U/t, U/s ^$> 1, where J\ = {t/s) 2 
and a — /3 2 , respectively. 



Let us discuss the 1/3 plateau problem by use of the effective Hamiltonian (TT9l) 
supposing J r ^> J\. In the 1/3 plateau case we can fix S z = +1/2 and retain the degree 
of freedom with respect to er, resulting in 



U 



(1/3) 
off 



tE( 






-(!-«)$>; 



(33) 



In other words, we take only two states ?/' ± (l/2, +1/2) in (fl8|) into consideration. We 



note that "Hg ff describes M = M s /3 states only, because both of -0^ =l= ' ) (1/2, +1/2) have 
Stot = +1/2 corresponding to M — M s /3. 

1/2 XY chain in transverse magnetic field, which 
1/2 XXZ chain in transverse magnetic field of which 



7i c g is nothing but the S 



is a special case of the S 
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ground state was investigated by Dmitriev et al [STMJT] . and Capraro and Gros 
Their results can read as follows in our cases. Both of (ofcj) and (— l) l ~^ (afa^) have 
the long-range order for the weak "magnetic field" case (i. e. < |1 — a\ <C Ji), 
while only (<7?er|) has the long-range order for the strong "magnetic field" case (i. e. 
|1 — a| ^> Ji). The boundary between above two cases will be approximately given 
by |1 — a\ ^ J\. In the latter case the state of the unit triangle is essentially either 
"00 (1/2, +1/2) or TpQ (1/2, +1/2) according as a < 1 or a > 1. In the former case, on 
the other hand, there exists the long-range Neel order of a spins along the x direction 
associated with the spontaneous breaking of the translational symmetry (SBTS) and 
the two-fold degeneracy of the ground state. In this case the ground state energy (the 
energy of the M s /3 state in the original spin language) will be considerably lowered 
by the interaction effects, which leads to the remarkable increase of the width of M s /3 
plateau. Further the SBTS of the a system results in the SBTS and the spontaneous 
breaking of the inversion symmetry of (S?) in the original spin representation. If we 
suppose the complete Neel order of o along the x direction (<rfj) = +1/2, (cr| J+1 ) = —1/2 
for simplicity (this situation will be a good approximation of the ground state in case 
of extremely weak "magnetic field" < |1 — a\ <C J\), the expectation values of S z are 

(Sf, 2j ) = (S z 3t2j+1 ) = 0.455 

(^, 2j ) = (^ + i) = - 167 ( 34 ) 

(Sfa) = (S Z 2W ) = -0.122 
which shows the spontaneous breaking of the translational symmetry (j =>- j + 1) and 
inversion symmetry (Sij <^> ^j)- 



(a) 1.0 




1.6 0.8 1.0 1.2 1.4 
a 




Figure 7. (a) (color online) a dependence of the M s /3 magnetization plateau width 
normalized by the saturation field H s for J\ — 0.1, 0.3, 0.5 cases obtained by DMRG 
with L = 128. The anomalous increase regions of A are shown by ellipses with red 
broken lines, (b) Behavior of the spin gap of the bond-alternating S = 1/2 Heisenberg 
chain described by the Hamiltonian (l36l) . 



Figure [7](a) shows the a dependence of the M s /3 magnetization plateau width A, 
from which we can see the anomalous increase of the A near a = 1 at least for J\ = 0.1 
and 0.3 cases as predicted by the above theoretical consideration. In usual cases the 
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plateau width (or the spin gap) is remarkably decreases near the mechanism-changing 
point and often becomes zero at that point, as shown in figure [3(b). Thus, we find 
a new and exotic behavior of the plateau width, which is completely opposite to that 
of usual cases [63]. In the J x = 0.1 case, for instance, the plateau is realized by the 
-0q (1/2, +1/2) and ipQ (1/2, +1/2) mechanisms for a < 0.9 and 1.1 < a, respectively, 
while by the cooperative effects of both mechanisms for 0.9 < a < 1.1. We note that 
the excitation gap (i. e., magnetization plateau width) A from the M = M s /3 state to 
M = M s /3 ± 1 states is always finite when a is varied. When J\ = 0.1 and a = 0.95, 
the expectation values (S z ) obtained by DMRG are 

(Si M ) = (Sfa +1 ) = 0.46 

(Sfo) = (Sl,2 j+ i) = 0-08 (35) 

(Sfa) = (S z 2,2 j+ i) = "0.04 
which well agree with the theoretically predicted values in (134]) . Details will be published 
elsewhere [58] . 

Let us explain why the new and exotic behavior of the plateau width is realized 
in our model, by comparing our model with the S = 1/2 bond-alternating Heisenberg 
chain described by 

H b _ a = J^{1 + (-1)^)5, • S j+1 (36) 

j 
where 8 is the bond alternation parameter. Any small amount of 8 produces the dimer 
spin gap as A oc \8\ 2 ' 3 [64ti66] (with leading logarithmic correction [67]). Thus the spin 
gap A of (136]) behaves as the sketch in figure [3(b) near 8 = 0. When 8 > (8 < 0), the 
spins S 2n and S 2n +i (>S^n-i an d S^n) effectively form a singlet dimer pair, which brings 
about the dimer spin gap. These two mechanism, 8 > and 8 < mechanisms, are 
completely competing and the reconstruction of the unit cell occurs at the mechanism 
changing point 8 = 0. At 8 = the system is of the Tomonaga-Luttinger state which 
is characterized by the gapless excitation. In our case, on the other hand, the unit 
cell is always the isosceles triangle for both of a < 1 and a > 1 cases. Then the 
reconstruction of the unit cell never occur in our case, which enables the cooperation 
of two plateau-formation mechanisms, ip^~\l/2, +1/2) mechanism and ip^ (1/2, +1/2) 
mechanism. 

The phase boundary between the traditional plateau phases and the new plateau 
phase where the plateau width remarkably increases can be determined by the 
phenomenological renormalization equation for the scaled gap LAo obtained by the 
numerically exact diagonalization method 

(L + 2)A 0iL+2 (« c , L ) = LA ,lK l ) (37) 

where A is the excitation gap within the M s /3 space and a c ^ is the size-dependent 
fixed point of a. The plateau-plateauless phase boundary can be estimated also by the 
phenomenological renormalization equation where Ao is replaced by A (plateau width) 
as in case of the spin gap at M = 0. Figure [HJ shows our preliminary results on the phase 
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Figure 8. Phase diagram of the M s /3 magnetization plateau. The plateau 1, 2 
and 3 correspond to the -0 (1/2, +1/2), the V> (1/2, +1/2), and the new plateau 
phases, respectively. Dashed lines are the second-order boundaries, while solid curves 
are expected to be of the Berezinskii-Kosterlitz-Thouless (BKT) transition. 



diagram at M s /3. This phase diagram qualitatively agrees with figure El However more 
detailed analyses will be necessary to draw a precise phase diagram. 

6. Twisted tube 

Although we have discussed theoretical aspects of the S = 1/2 isosceles quantum 
spin tube, a compound corresponding to our model has not been found yet. In this 
section, we review the twisted spin tube, which is of particular importance for actual 
experiments; Recently, the interesting compound [(CuCl2tachH) 3 Cl]Cl2 was actually 
synthesized as an assembly of the triangular cluster conformed by Cu ++ ions and 
magnetization measurements were performed |2]. As is in figure El the tube structure 
of this compound is based on the alternatingly-aligned triangles, where the exchange 
coupling in the unit triangle is denoted as J r and the intra-triangle coupling is written 
as 3\. The low-energy properties of the quantum spin tube were theoretically studied 
by DMRG [HJ[T5]. Also, the quantum spin tube with the easy-plane anisotropy was 
investigated in the context of the triangular lattice spin system and then interesting 
field induced phase transitions were observed in the magnetization curve |68j . 

At the early stage of study, [(CuCl 2 tachH) 3 Cl]Cl 2 might have a small spin gap, 
based on the susceptibility measurement above 2K. However, the succeeding studies 
have clarified that the observed spin-gap-like behavior around 2K crossovers to the 
gapless behavior at very low temperature region (T < 0.5K) [69], which is consistent 
with the theoretical results. This suggests that a kind of dual structure associated with 
the quantum phase transition of the first order may be captured by the experimental 
measurements. Below, we will review the quantum phase transition peculiar to the 
twisted quantum spin tube. 
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Jl 




Figure 9. Twisted quantum spin tube. J r indicates the exchange coupling in the unit 
triangle and J\ does inter triangle coupling. The expansion of the twisted spin tube 
has the triangular lattice structure. 



The model Hamiltonian of the twisted spin tube is written as, 

3 L 2 L 

t-l = Ji 2_^ 2_^ [Si,j ■ Sz,j+i + Si,j ■ Si+ij+i] + j r y j y j s itj ■ s i+ ij (38) 

i=l j=l i=l j=l 

with the periodic boundary condition in the z-direction. A naive consideration about the 
lattice structure of this twisted tube suggests that a quantum phase transition between 
the following two phases may occur; For J r ^> J\, the system is described by the almost 
decoupled triangles, where the chirality degree of freedom defined on the unit triangle 
plays crucial role, while for J r <C J\, the lattice structure is basically described as a 
rhombus lattice. Note that, if J = 0, the rhombus lattice has no frustration. Thus 
we can expect a quantum phase transition between decoupled-triangles and rhombus 
lattice phases. 

In the decoupled triangle limit (J r 3> Ji), the degenerate perturbation with respect 
to Ji leads us to the effective spin-chirality model, 

2 7 

U ^ = X E T i ■ T ^ 1 + 2(exp(i7r/3)r/r- +1 + (exp(-in/3)r7 T + +1 )} (39) 

3 

where T represents the S = 1/2 spin operator and r indicates an effective spin matrix 
representing the chirality degrees of freedom as in (j2"2"|) . This effective Hamiltonian is 
very similar to f[2"2"j) ; The phase factor e 47r//3 originates from the 7r/3 rotation of the 
lattice along the tube direction and it can be removed by the gauge transformation. 
On the same line of the argument as ( J22l) . the effective Hamiltonian ( )39l has a spin 
gap [ini[IIl[13]. Indeed, a detailed DMRG computation of the full Hamiltonian (1381) 
actually confirms that the spin gap exists for Ji/J r < 1.21 [HI [15]. F° r Ji/Jr > 1-21, 
on the other hand, the finite size scaling analysis of the spin gap basically indicates 
the gapless ground state [15]. Although the usual singlet-triplet spin gap for the open 
boundary condition captures a boundary excitation in 1.21 < Ji/J r < 1.5, such an 
anomalous behavior due to the boundary can be settled down by using the single-spin 
termination of the tube [8]. 

In order to resolve nature of the quantum phase transition of the twisted tube, an 
important quantity is the total-S" on the unit triangle [BJ. For decoupled triangle, 
the ground state is in doubly degenerating doublets of the unit triangle, while for 
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Figure 10. (color online) The expectation value of the projection operator (-P1/2) for 
the extended quantum spin tube. The solid line represents 7 = and the solid circle 
indicates the parameters corresponding to [(CuCl2tachH)3Cl]Cl2. 



J r = 0, the ground state basically belongs to S — 3/2 sector. Here, we define 
the projection operator of the total-S" of the unit triangle to the doublet sector as 
Pi/ 2 (= \ — \{S\ ■ S2 + S2 ■ S3 + S3 • Si)) for each unit triangle. As discussed in [H], 
then, this projection operator Py 2 is a good order parameter of the twisted tube; The 
detailed computation of the expectation value (-P1/2) shows discontinuity at the critical 
value (Ji/J r )c — 1-22 not at J\/J r ~ 1.47, which is consistent with the spin gap result. 
In addition to this, we should remark that an extended spin tube with the diagonal 
interaction shows the exact first order phase transition, for which total-S 1 on the unit 
triangles independently conserves [70J. We have thus investigated how to connect the 
twisted tube to the diagonal interaction model. We add 

3 L 

KjJi = 1-h ^2 5Z S U ' S i-U+i ( 40 ) 

i=i 3=1 

to the Hamiltonian (1381) . and then varies — 1 < 7 < 1. Note that 7 = corresponds to 

the original twisted tube and 7 = 1 does to the diagonal spin tube, which shows the 

first order transition at J\J J t ~ 0.63 [70] . 

In figuredUl we show (-P1/2) at the center of the extended tube calculated by DMRG 

for L = 36 with J r = 1. In the figure, (-P1/2) jumps exactly from zero to unity at 

J x j J t ~ 0.63 along the line of 7 = 1, where the total-S" of the unit triangle is the exact 

symmetry. Then, we can see that the "gap" of the (-P1/2) is adiabatically continued 

to 7 = 0. This implies that the transition at {J\/ J T ) C — 1-21 on the line of 7 = 

corresponding the twisted tube fl38l) is of first order. As J\ increases in 7 < 0, this gap 

reduces, but we can not confirm whether the end point of the first order transition exist 

or not within the present calculations for L = 36. However, we think that the overview 

of the first order quantum phase transition is clarified. 
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In connection with the experiments, an interesting point is that the parameter 
Ji/ J r ~ 2.16 for [(CuCl 2 tachH) 3 Cl]Cl 2 is located in the gapless S = 3/2 sector, but it is 
surrounded by the "gap" of the first-order transition, as in figure fTUl This implies that, 
although the ground state itself is gapless, a certain density of states originating from the 
S = 1/2 sector can be expected above the ground state. Thus we can expect that such 
dual structure of S — 1/2 and 3/2 sectors affects the experimentally observed quantities 
of the twisted tube at a finite temperature. Indeed, a recent specific heat result illustrates 
that the linear temperature dependence at very low temperature(T < 0.5K), which can 
be described as an effective S = 3/2 chain, crossovers to spin-gap-like increasing around 
T ~ 2K [69]. Also for the low-field magnetization curve of J\/ J x > 1.21, we can see 
that, after linear increase of the magnetization at very low field region, the slope of the 
magnetization curve rapidly increases, as if it had a spin gap [15] . 

7. Future prospects 

As a future prospect, it would be interesting to consider the carrier-doped spin nanotube 
like the high-Tc cuprates, where the system is effectively described by the Hubbard 
model near the half-filling case. A mechanism of the superconductivity based on the 
spin gap had been proposed for the carrier-doped two-leg spin ladder [71] and actually 
a pressure- induced superconductivity was observed on the spin ladder cuprate [72J. 
Motivated by the discovery, the spin gap mediated superconductivity was theoretically 
proposed even for the three-leg spin ladder, using the quantum Monte Carlo simulation 
[?3fT4"] . If the carrier-doped three-leg spin tube is realized, it would be a better candidate 
of the superconductor, rather than the three-leg ladder, because it has a spin gap. In 
addition the three-leg spin tube was also revealed to have an energy gap in the chirality 
degrees of freedom [13J. It suggests that a chirality induced superconductivity would be 
possibly realized as a new mechanism of superconductivity in the near future. 
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